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ABSTRACT 

The stellar haloes of galaxies can currently be studied either through observations of 
resolved halo stars or through surface photometry. Curiously, the two methods appear 
to give conflicting results, as a number of surface photometry measurements have 
revealed integrated colours that are too red to be reconciled with the halo proper- 
ties inferred from the study of resolved stars. Several explanations for this anomaly 
have been proposed - including dust photoluminescence, extinction of extragalactic 
background light and a bottom-heavy stellar initial mass function. A decisive test is, 
however, still lacking. Here, we explain how observations of the halo of a nearby galaxy, 
involving a combination of both surface photometry and bright star counts, can be 
used to distinguish between the proposed explanations. We derive the observational 
requirements for this endeavour and find that star counts in filters VI and surface 
photometry in filters VIJ appears to be the optimal strategy. Since the required halo 
star counts are already available for many nearby galaxies, the most challenging part 
of this test is likely to be the optical surface photometry, which requires several nights 
of exposure time on a 4-8 m telescope, and the near-IR surface photometry, which is 
most readily carried out using the upcoming James Webb Space Telescope. 
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1 INTRODUCTION 

The stellar haloes of galaxies contain the fossil record of 
galaxy assembly, and observations of resolved halo stars 
in th e M ilky Way, Andromeda and other nearby galaxies 
(e.g. iMou hcinc ct al. 2005; Ibata ot al. 2007; Hclmi ct al. 



20081: iReikuba. Mouchine fc Ibatall2009i : iDurreU et al.ll2O10l : 
Tanaka et aLlbOHl ) have provided a wealth of information 
on such structures. However, stellar haloes can also be stud- 
ied through surface photometry, i.e. observations of the in- 
tegrated light from large numbers of unresolved stars within 
each system. 

Both techniques have their strong sides as well as short- 
comings. Star counts can trace structures much further out 
into the halo, but this method is currently limited to the 
brightest halo stars (with the Milky Way halo as a notable 
exception) and to systems within ~ 10 Mpc. Surface pho- 
tometry can be applied for more distant galaxies, is sensitive 
to less luminous halo stars and diffuse light from the inter- 
stellar medium, but works only for the relatively bright inner 
regions of halos. 

At the current time, the results from these two tech- 
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niques are difficult to completely reconcile. A number of at- 
tempts to study the haloes of edge-on disk galaxies through 
optical and near-infrared (NIR) surface photometry have re- 
vealed halo colours that are too red to be consistent with 
the halo po pulations captured through the study of resolved 
stars (e.g. 



Rudy et al 



Sackett et al 
19971: 



1994'; 'Lequeux et"al] 19951. 



James fc Casali 1998: 



Zib etti et al.ll2004l : Tzibetti fc Fergusonj 



Lequeux et al 



1996 



1998 



^ | |2004 iBergvaU et al.l 

,20ld r The most recent observations indicate that this 
red excess turns up at extremely fa int surface bright- 
ness levels (ni « 27-29 mag arcsec"^; IZibetti et al. I l2004l : 
iBergvaU et all l2010l i and at wavelengt hs around the ii- 
band (» 7600 A). O lder measurements l|Rudy et al.lll997l : 
I James fc Casllil[l998l ) suggest that it may also continue into 
the near-infrared (12000-22000 A). On the other hand, sur- 
face photometry carried out at wavelengths sh ortward of the 
ii-band seem to turn up nothing anomalous ( Zibetti et al.l 
|2004 [ Jablonka et al. I I2OIOI : iBergvaiT et al. I l2010l v 



It has been suggested that these 'red halos' could be 
due to mundane obse rvational effects like instrumental scat- 
tering (|de Jondl2008h or flux residuals left over from an in- 
complete sky subtraction l|Zheng et al. I ll999l : I Jablonka et al.l 
|2010D . In some of the claimed detection s, this may well be 
the case. However, IBergvaU et"all (|201Gl i correct for instru- 
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mental scattering (point spread function) effects and are still 
left with a pronounced red excess in their study of halos 
around stacked, edge-on low surface brightness disks. They 
also show that while the halos of these disks have extremely 
red colours, the disks themselves do not display any anoma- 
lous colours at the corresponding surface brightness levels. 
This strongly argues against a sky flux residual as the main 
culprit, since this would result in strange colours in both 
components. 

The red halo phenomenon is not necessarily confined 
to disk galaxies only. Near-infrared excesses have been 
rep orted in the faint out s kirts of blue comp act galax- 
ies llBergvall fc Osthnll2002l : iBergvaU et "all 120051), alt hough 
IZackrisson et al.l l\20m ) and iMicheva et al.l (|2010l ') have 
questioned the need for extreme stellar population proper- 
ties to explain these. There are also older reports of ex- 
tremely red op tical colours in the outskirts of cluster ellip - 
tical galaxies l|Maccagni et al.|[l98l : iMolinari et al.lll994 ). 
However, no similar anomali e s wer e seen in the recent 
study by iTal fc van DokkumI (|2010l ) of the outskirts of 
42000 stacked early-type galaxies (mostly located in galaxy 
groups) . 

If one assumes that the red excess in the halos of edge- 
on disks is not due to instrumental scattering or sky flux 
residuals, this still leaves room for a number of possible ex- 
planations - including a bottom-heavy stellar initial mass 
function, diffuse light from the interstellar medium and ex- 
tinct i on of extragalactic backgr ound light (I Zackrisson et all 
l2006l . l2009l : [BergvaU et alll2010l ). These solutions can all give 
rise to similar halo colours. Given the extremely challenging 
task of improving current surface photometry measurements 
with existing telescopes, it has so far remained unclear if the 
degeneracy between these really can be broken. In this pa- 
per, we describe a method to discriminate between these 
scenarios using a combination of star counts and surface 
photometry. 

In Section. [2j we review a number of explanations for 
the red excess that have previously been demonstrated to be 
broadly consistent with the available data. The principle of 
our proposed test is described in Section |3] and the observa- 
tional requirements derived in Section|3] Section[S]discusses 
a number of potential target galaxies that meet the distance 
requirements of the test. We also provide a worked example 
for the required exposure times in the case of the edge-on 
disk galaxy NGC 4565, which is deemed to be one of the 
most suitable candidates. Our findings are summarized in 
Section |B1 



2 POSSIBLE EXPLANATIONS FOR THE RED 
EXCESS 

2.1 A bottom- heavy stellar initial mass function 

It has long been suspected that the red excess could be a sig- 
nature of a stellar halo wit h an extremely bottorn-heavy stel 



1997 




2ooe 





Zepf et al. 2000: Z ibctti ct al. 2004; Zackris son et al 
Bergvall et al.1 l20"ld ). i.e. a stellar halo with an ab- 



an explanation for the dark matter of galaxies. However, 
it is now well established that the dark ma tter is predom- 
inantly non-baryonic on cosmic scales (e.g. iKomatsu et al.l 
[2011 ). Surveys for microlensing events from Massive Astro- 
physical Compact Halo Objects (MACHOs) have also ruled 
out a dominant contribution to the dark h alo of the Milky 
Way from objects in this mass range (e.g. iTisserand et al.l 
2007*). On the other hand, up to « 2/3 of the cosmic 
baryons in the low-redshift Universe a re still unaccounted 
for (e.g. IProchaska fc TumhnsonI |2009| ). and the missing- 
baryon problem may be even more severe o n scales of in- 
dividual galaxies (e.g. iMcGaugh et al.ll2O10l ). There is also 
a long-standing MACHO puzzle, in the sense that some 
halo microlensing surveys have turned up more events than 
can easily be attributed to known stellar populations (e.g. 
Alcock et al."2000'; 'Calchi Novati et aLllioOsI : llngrosso etal] 
_2007; Riffcscr, Scitz & Bender 2003). Hence, it could be ar- 
gued that there may be low-mass stars hiding in the halo 
that contribute to the missing-baryon problem and/or to the 
MACHO microlensing event rates, even though they have 
little bearing on the dark matter. 

The stellar initial mass function (IMF) has long 
been assumed to be universal, but recent observa- 
tional studies based on unresolved stellar popula- 
tions have given some support to the n otion that it 
may vary as a functio n of cosmic time (van Dokkunil 
2008). environment (iHoversten fc Glazebrookl l2008l : 



Meurer et al] l2009l: iLee et al.l |2009| ) and star formation 



rate I Gunawardhana et al.l I2OII ). The strength of cer- 



normally high fraction of low-mass stars (with masses < 
0.5 Mq). Due to the high mass-to-light ratios of such popu- 
lations, solutions of this type were once thought to provide 



tain absorption features in giant elliptical galaxies have 
been interpreted as evidenc e of a bottom-heavy IMF 
(Ivan Dokkum fc Conrovl [2OI0I ), whereas the same absorp- 
tion lines are cons istent with a standard IMF in the case of 
globular clusters (|van Dokkum fc Conrovl I2OIII ). Attempts 
to constrain the IMF in the field popula tion of the Large 
Magellanic Cloud using star counts (e.g. ICouliermis et ahl 
l2008i) have also yielded results seemingly inconsisten t with 
the standard IMF (e.g. lKroupall200ll : IChabrieij|2003l ). 

Surveys for low-mass stars (subdwarfs) in the Milky 
Way halo have produced n o evidence for a bottom-heavy 
IMF (e.g. lCould et al.lll998l '). but current subdwarf searches 
are limited to distances out to « 40 kpc from the Sun. If 
there were a second halo population with a bottom-heavy 
IMF and a density profile different from that of the hith- 
erto known halo, this second population could have evaded 
detection in the subdwarf searches. Such a halo population 
could on the other hand turn up in surface photometry mea- 
surements of external galaxies (and thus explain the red ex- 
cess), due to projected contributions from stars further away 
from the halo centre . Scenarios of this type were explored by 
IZackrisson fc Flvnnl l|2008l ). who used the existing subdwarf 
data to rule out all secondary halo populations that would 
be able to contribute significantly to the missing baryons. 
Moreover, a halo like this would only remain viable if the 
surface brightness of the red halo reported bv lZibetti et al.l 
l|2004h for stacked edge-on disks had been significantly over- 
estimated. This could for instance happen because of instru- 
mental scatt ering, which Zibetti e t al. d id not properly cor- 
rect for. The IZackrisson fc Flvnnl (I2OO8I ) constraints, which 
are much stronger than those currently produced by mi- 
crolensing searches, could also be sidestepped if the low- 
mass stars in the second halo were all locked up in star 
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clus ters. This, on the oth er hand, would cause tension with 
the iBergvaU et all l|2010l ) stacked halo observations, which 
favour small pixel-to-pixel variations in halo colours (i.e. a 
rather diffuse and not too highly clustered halo). 

In summary, there are no observations that firmly rule 
out a halo component with a bottom-heavy IMF, but ex- 
isting microlensing data, subdwarf star counts and surface 
photometry measurements do impose strong constraints on 
the properties of such a population. Current data place the 
mass that can be locked up in such a population at a level 
that is too low (<C 10^" Mq) to have any significant bearing 
on the missing baryon problem (~ 10^° Mq), and definitely 
not on the dark matter prob lem lO'^^ Mp)) on the scale 
of Milky Way-sized galaxies (jZackrisson fc FlvnulbOOSl ). 



2.2 Diffuse light from the interstellar medium 

Another possibility is that the red excess does not origi- 
nate from the stars in the halo, but is instead caused by 
some process in the interstellar medium. Diffuse, low-density 
halo gas could be kept ionized by far-ultraviolet photons 
leaking from the disk or by shocks and would contribute 
emission lines and a nebular continuum to the integrated 
halo spectrum. This could indeed give rise to strange halo 
colours, but current predictions for superpositions of old, 
low-metallicity stellar populations and ionized halo gas are 
still unable to match the observed colou rs of red halos a t 
low redshifts (Zackrisson et al. 2006; Bc rgvall et al.|[201Cf ). 
Due to the high emission-line equivalent widths of a purely 
nebular spectrum, the optical colours of ionized halo gas 
ev olve dramatically wi t h red shift, and the colours reported 
bv lZibetti fc Ferguson' (12004 ) for a halo around an edge-on 
disk galaxy at redshift z = 0.322 can in principle be ex- 
plain ed by a mixture of ion ized gas and a normal stellar 
halo l|Zackrisson et al1l2006l ). However, the halo colours of 
this particular target may also be s everely affected by in- 
strumental scattering (|de Jong||2008l ). 

If the red excess originates in the interstellar medium, 
dust emission seems a more plausible radiation mechanism. 
Mid-infrared dust emission has bee n detected in the halos 
of several edge-on disk galaxies (e.g. Ilrwin fc MaddenI l2006l : 
iBurgdorf et~al]|2009l : iKaneda et al.ll2009^■ but to account for 
th e red r — i halo colou rs reported bv IZibetti et al.l (|2004l ) 
or iBergvall et aD (|2010l ). the dust emission would need to 
extend into the optical i-band. In principle, this could be 
achieved by the dust photoluminescence phenomenon known 
as extended red emission (ERE). While the dust compo- 
nent responsible for the ERE has yet to be identified, ERE 
has been reported in a wide range of environments includ- 
ing refiection nebulae, HII regions, the halo region of the 
starburst galaxy M82 and the high-latitude diffuse interstel- 
lar medium and c irrus in the Milky Way (for a review, see 
I Witt fc Viihll2004h . ERE manifests itself as a wide (~ 1000 
A) emission feature in the wavelength range from 6100-9500 
A, where the peak wavel ength appears t o vary from object 
to object. As argued bv iBergvall et al.l ([2010), ERE can- 
not be excluded as the mechanism responsible for the red 
r — i halo colours, but an omalous colours involving filters 
at longer wavelengths (e.g. iRudv et al.lll997l ) would then re- 
quire a second explanation. 



2.3 Extinction of extragalactic background light 

Anomalously red colours (as measured by colour indices like 
V — I,r — ioxV — K) do not necessarily imply the pres- 
ence of a flux excess in the longer wavelength filter (e.g. 
the i band in the r ~ i colour index), since red colours can 
also be produced by a flux deficit in the shorter wavelength 
band (e.g. the r band in r — i). This could most readily be 
caused by dust reddeni ng of the halo stars, b ut as shown by 
[Zackrisson et all (|2006l ) and IBergvall et al ] (2010'), the ob- 
served red halo colours do not lie along standard dust red- 
dening vectors. Moreover, the reddening needed to explain 
the observed r — i halo colours would require A{g) « 1-3 
magnitudes of p-band extinction, which seems completely 
implausible at projected distances of ~ 5-10 kpc into the 
halo. 

However, as de monstrated by IZackrisson et al.1 (Hooi) 
and IBergvall et al] (|2010l ). there is one mechanism that 
could account for the observed halo colours while requir- 
ing no more than A{g) ~ 0.03-0.2 mags of optical ex- 
tinction, which would be consistent with the halo opaci- 
ties inferred from the reddening of background light source s 
(IZaritskvlll994l : iMenard et al]l2010l : iMcGee fc Baloghll2010l l . 
This mechanism is based on the extinction of extragalactic 
background light (EBL) , which gives rise to a systematic sky 
subtraction problem for surface photometry measurements 
at faint isophotes. While most of the sky flux originates from 
regions between the telescope and the haloes targeted by 
surface photometry observations, a small fraction - the EBL 
- comes from behind. Provided that the ex isting direct mea - 
surements of the optical EBL are correct (|Bernsteinll2007l l. 
most of this light appears to be diffuse (i.e. unresolved with 
all existing instruments). Unlike the other components of the 
night sky flux, the EBL can be subject to extinction by dust 
present in the target halos. This invalidates an assumption 
adopted in all current surface photometry measurements, 
namely that the sky flux and the flux from the target ob- 
jects are unrelated. This extinction of the EBL may lead to 
a slight decrease in the overall sky flux across the target ha- 
los, which is neglected when the sky flux level is estimated 
well away from the target objects. The net effect is a sys- 
tematic, wavelength-dependent oversubtraction of the sky, 
with spurious colour gradi ents in the outskirts of extended 
targets as the likely result (|Zackrisson et al1l2009l ). 

If the red halo colours can be demonstrated to be due 
to EBL extinction, this could provide interesting constraints 
on the EBL itself. For halo opacities in the A{g) ~ 0.01- 
0.1 range (in line with current measurements), the surface 
brightness of the optical EBL would have to b e at a level 
similar to that suggested by the lBernsteinI (|2007h EBL mea- 
surements. This is considerably higher t han what curren tly 
resolved galaxies can account for (e.g. iMad au fc Pozzettil 
I2OOOI I. Hence, this would imply that the origin of the op- 
tical EBL remains unknown, with potentially far-reaching 
implications for cosmology. 



3 A DECISIVE TEST 

3.1 Relating surface photometry to star counts 

In recent years, most attempts to trace the surface bright- 
ness profiles of stellar haloes have been based on direct star 
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Figure 1. Conversion between total surface brightness /^total (as measured through surface photometry) and the surface brightness due 
to resolved stars ^counts for an instantaneous-burst (i.e. single-age) stellar population, a) The /-band conversion factor /^total ~ /^counts 
as a function of age for metallicities Z = 0.0001 (thin line), Z = 0. 001 (med i um th ick line) and Z = 0.008 (thick line), in the case where 
stars are resolved 2.0 mag below the TRGB. Solid lines are for a ISalpetej ([1955') IMF between 0.08 and 120 Mq, whereas the dashed 
lines corresponds to the iKroupal ll200ll ) IMF corrected for binary stars. The difference in /^total ~ /^counts between these two IMFs is 
insignificant, but variations with age and metallicity variations are not. Still, the variations with age are small for old ^ 6 Gyr stellar 
populations, and the metallicity can typically be estimated from the colours of the TRGB, which means that uncertainties due to the 
metallicity can be corrected for. b) The difference between the integrated V — I colours of all resolved stars and the integrated V — I 
colours of the entire population, {V — /)tot — {V — /)counts, in the case where stars are resolved 2.0 mag below the TRGB. The lines have 
the same meaning as in panel a. The variations with age are relatively small (~ 0.2 mag), but the variations with metallicity remain 
substantial (si 0.5 mag), c) Same as b, but for counts reaching 4.0 mag below the TRGB. The uncertainty in /^total — /^counts due to both 
age and metallicity variations is significantly reduced, d) Same as b, but for counts reaching 4.0 mag below the TRGB. The sensitivity 
of (V — /)tot — {V — /)counts to age variarions is very small, but metallicity variations can still produce significant colour shifts 0.2 
mag). 



counts, i.e. by counting the number of bright halo stars (red 
giant branch stars, asymptotic giant branch stars and upper 
main sequence stars) within a certain area. 

The bright stars only account for a fraction of the over- 
all halo hght, which means that the star counts must be 
calibrated onto a surface brightness scale by making assump- 
tion about the stars below the detection threshold. Usually, 
this is achieved by simply shifting the star counts profile 
to match actu al surface brightness measurements a t bright 
isophotes (e.g. llrwin et al.ll2005l : [Barker et al.ll2011^ . For an 
old stellar population with a reasona bly normal IMF (e.g. 
ISalpeteil[l955l : lKroupall200ll : IChabrier 2003), this procedure 
is fairly robust, as demonstrated in Fig. [1] where we present 



the conversion factor as a function of age, metallicity, IMF 
and depth of the star counts. These conversi on factors are 
derive d by integrating fluxes of stars along the lMarigo et al] 
(|2008l ) stellar isochrones (with suitable IMF weights ap- 
plied), and comparing the surface brightness /^totai produced 
when the entire isochrone is used (the case relevant for sur- 
face photometry) to the surface brightness /icounts produced 
when only stars above a certain flux threshold are included 
(the case relevant for star counts). 

To illustrate how this conversion factor depends on the 
depth of the star counts, two cases are explored - star counts 
reaching 2.0 mag below the tip of the red giant branch 
(TRGB) in the /-band (this corresponds to a detection 
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Figure 2. Diagram illustrating how the possible explanations for red halo colours can be tested using a combination of star counts and 
surface photometry, a) The halo surface brightness profile in V (blue) and / (red) as inferred from surface photometry (lines) and star 
counts (filled circles) in the case where the red halo colours are caused by excess flux in both the V and / bands, as would be expected 
from a stellar population with a bottom-heavy IMF. b) Same in the case where the red halo colours is caused by dust photoluminescence 
(ERE), resulting in a flux excess limited to the / band, c) Same in the case where the red halo colours are caused by a larger flux deficit 
in the V-band than in /, due to EBL extinction, d) The colour profile resulting from either of the situations depicted in a) and b) and 
c). See main text for additional details. 



threshold of Mi « -2) at an age of 10 Gyr in a Z = 0.001 
population (Fig. [T^ & b), and star counts reaching 4.0 mag 
below the TRGB (detection threshold M/ ~ 0) for the same 
age and metallicity (Fig. [TJ; & d) . For simplicity, we assume 
that the star counts are sensitive to all stars brighter than 
the cutoflF, and do not count stars in diflterent evolved evolu- 
tionary stages separately. The TRGB provides a convenient 
reference point as its /-band magnitude of A/j ~ —4 is nearly 
independent of age and metallicity for any population older 
than 3 Gyr and more metal-poor than 0.5 x Zq, hence its 
use as a distance indicator of halo-like stellar populations. 



In summary, even if the s tellar halos of galaxies contain 
substantial age grad ients (e.g.|Stinson et al.i .2009) or metal- 
licity gradients (e.g. ICarollo et al.l 20071 ). the /itotai — Mcountss 
uncertainties are small, provided that the integrated light 
comes from stars only, that the metallicity can be obser- 
vationally assessed and that the IMF is of standard form. 
However, the reports of anomalously red halo colours based 
on surface photometry indicate that one (or several) of these 
tenets may be violated. 



While /itotai — /^counts shows substantial variations (~ 1 
mag) as function of both age and metallicity in the case 
of the shallower star counts (Fig. [T^), the variations with 
age are small (~ 0.2 mag) in the age range relevant for old 
stellar haloes (> 6 Gyr) if one assumes that the metallic- 
ity is known. The /xtotai — ^'counts variations with age are 
moreover about the same in V and /, so that the diflterence 
between the integrated colours of all resolved stars and the 
integrated colours of the entire population become no more 
than (y-7)tot-(l^-/)counts « 0.1 at ages > 6 Gyr (Fig.fljj). 
While this leaves room for substantial uncertainties in both 
/itotai - /^counts and (V - /)tot — {V - /)t:ounts due to ra- 
dial metallicity variations, such metallicity gradients would 
result in radially varying colours for the TRGB, and can 
therefore be observationally tested for. For the deeper star 
counts (Fig. [l}; & d), both uncertainties are substantially 
smaller. 



3.2 Profile offsets 

The solutions to the red halo puzzle reviewed in Section [5] 
can be separated into two broad categories, depending on 
how the surface brightness profile inferred from surface pho- 
tometry would diflter from that inferred from star counts. 
Bottom-heavy IMFs and dust photoluminescence (ERE) 
would both decrease /itotai — Mcounts (i.e. make the surface 
photometry profile brighter than the star counts profile), 
whereas EBL extinction would tend drive this quantity in 
the opposite direction. The principle is schematically illus- 
trated in Fig. [21 where we show how the halo surface bright- 
ness profiles based on star counts and surface photometry 
would deviate from each other in these three cases. Here, 
we have assumed that the star counts profiles have already 
been shifted to the level of the surface brightness profiles 
at bright isophotes (where a standard stellar population is 



6 Zackrisson et al. 



0.5 



t 



Bottom-heavy IMF, V-band 
Bottom-heavy IMF, l-band 
Extended red emission, V-band 
Extended red emission, l-band 
EBL extinction, V-band 
EBL extinction, l-band 




4 6 
R (kpc) 



10 



Figure 3. The offset Afi between the total surface brightness (as 
measured through surface photometry) and that inferred from 
star counts, in filters V (blue lines) and I (red lines) in the 
case where either a bottom-heavy stellar IMF (solid lines), ERE 
(dashed lines) or EBL extinction (dash-dotted lines) are respon- 
sible for pushing the integrated V — I colours to anomalously red 
values (Fig. [2ji). A positive value on Afj, implies that the sur- 
face photometry profile lies faintward of the star counts profile, 
whereas a negative value implies that surface photometry profile 
lies on the brightward side. The three scenarios result in qualita- 
tively different predictions in this diagram. 



assumed to dominate), so that the /itotai — ^counts offsets 
evident from Fig. [1] have already been corrected for. 

In the first case (Fig. [2^), we assume that there is a halo 
component with a normal V — I colour (similar to that of 
the stellar population hitherto detected in the Milky Way 
halo) that is slowly overtaken by a bottom-heavy IMF pop- 
ulation with a much redder colour, thereby driving the inte- 
grated light measurements to extreme V — I values at large 
projected radii. Since bright star counts fail to detect the 
increasing number of low-mass stars in the outer profile, the 
total surface brightness is shifted in the brightward direc- 
tion in both the V and the / bands, although more so in the 
/-band. 

In the second case (Fig. ^jp) , we assume that the same 
standard halo is affected by ERE in the /-band at large 
projected radii, thereby driving the integrated V — I to the 
same extreme values as in Fig. [5^. Whereas surface pho- 
tometry may agree with the star counts in the V-band, the 
/-band surface photometry data would indicate a higher sur- 
face brightness than that inferred from star counts. 

In the third case (Fig. [2};), the colours of the normal 
halo component are pushed towards extreme values due to 
a failure to account for EBL extinction in the surface pho- 
tometry measurements. Both the V and / surface photom- 
etry profiles would then indicate a surface brightness lower 
than that inferred by star counts. Here, the anomalously red 
integrated colours arise from a greater fiux deficit in the V 
band than in /. 

The model parameters have been chosen to produce the 
same integrated V — I profile (Fig. Oi) in all three scenarios, 
starting from V — I = 1.15 (typical of an old, metal-poor 
stellar population) at small radii and reaching V ~ I — 1.45 



(anomalously red) at the outermost radius plotted. The off- 
sets A/i between the surface photometry and star counts 
profiles in Fig. [5] have been plotted in Fig. O As seen, the 
directions and magnitudes of the V and / offsets are qualita- 
tively different in these three cases. At the outermost data 
point, where a colour of V — I 1.45 has been reached, 
the amplitude of the offset is 0-0.5 mag in the V-band and 
0.2-0.75 mag in the /-band. 

We argue, that by combining star counts and surface 
photometry for a single target in filters V and /, one should 
at least be able to distinguish between an excess of red (/- 
band) fiux (bottom-heavy IMF or ERE, as in Fig. [2^ and 
b) and a deficit of blue (l/-band) fiux (EBL extinction, as 
in Fig. [5^) as the reason for the red halo colours. While 
Fig. |3] suggests that V and / data should also be able to 
distinguish between a bottom-heavy IMF from ERE, based 
on the lack of offset in the l/-band in the latter case, this 
conclusion would not be as robust. While Fig. |3] certainly 
depicts a posstble signature of ERE, it is not necessarily 
generic, simply because the colours of ERE are not suf- 
ficiently well-constrained. The ERE emission peak, which 
we have assumed to affect the /-band only, tends to co- 
exist with a blue component of dust-scattered continuum 
radiation, whose relative contribution depends on the angle 
between the line of sight and the direction from which the 
ERE- fuelling ultraviolet photons are originating (in our case 
supposedly the disk). If this diffuse continuum is sufficiently 
strong, the offset Afi between the surface photometry and 
star counts profiles would become non-zero in the V band 
as well, thereby producing a situation similar to that of a 
bottom- heavy IMF (Fig. [gja,). 

However, if the outcome favours an excess of red fiux 
{A^v and A^i negative; |A/ii| > lA/iv'l in Fig. |3} over a 
lack of blue fiux (A/iv and A^i positive; |A/ii| < |A/iv| in 
Fig. [3} as the reason for anomalous halo colours, a bottom- 
heavy IMF can be distinguished from ERE by additional 
surface photometry data in a filter at longer wavelengths 
(J, H or K). In the case of a bottom-heavy IMF, the fiux 
excess should grow stronger towards longer wavelengths, but 
not in the case of ERE, since this emission peak is limited 
to the wavelength range 6100-9500 A. Due to observational 
limitations (see Section |4ll , the J band appears to be the 
most suited for this endeavour. 

There is admittedly one possibility that is not covered 
by the scenarios presented in Figs. [5] and O Let's assume 
that the seemingly abnormal halo colours are caused by a 
failure of the models to adequately reproduce some bright, 
red phase of stellar evolution. This omission would need to 
be present in all of the many models alread y applied to the 
red halo problem (e.g. iBergvall et al.|[2010l ), and since the 
halo colours are considerably redder than the observed to- 
tal colou rs of old systems lik e elliptical galaxies or globular 
clusters jZibetti et al.|[200^ it is perhaps more likely that 
there is a problem with the models at low to intermediate 
ages. If the integrated light of some region of the halo would 
be dominated by a relatively young, red component, the ac- 
tual value of /itotai — Mcounts could be closer to zero than 
predicted in Fig. [TJ since a larger fraction of the integrated 
light would come from resolved stars. Provided that this 
happens at a radius greater than that at which the /itotai 
and /icounts are normalized to the same value, this would 
cause the star counts (dots) to rise above the surface pho- 
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tometry profile (solid line) in Fig. (2] Unlike the case for 
EBL extinction (Fig. [2};) , the rise in the star counts would 
be more prominent in / than in V (A^y and A^i positive; 
|A/x/| > |A/iyl in Fig. [3J, because of the very red colours 
of these bright objects. If this is accompanied by colour ex- 
cess in the surface photometry measurements, this would 
be a tell-tale signature that the orgin of the abnormal halo 
colours should be sought among the resolved objects. How- 
ever, a situation of this type is in conflict with the small 
di fference between mean and median halo stacks measured 
bv lBergvall et al.l (|2010| ). which suggest that the sources re- 
sponsible for the red halo colours are smoothly distributed 
throughout the halo and not concentrated to a few bright 
pixels. A sudden excess of very red, bright objects at some 
radius should also be directly detectable in the halo colour- 
magnitude diagram. 

The surface brightness profiles depicted in Fig. [2] have 
been generated assuming that the surface brightness profile 
of the standard halo, as measured by star counts (corrected 
for the /itotai — Accounts offsets in Fig. [1} is described by a 
Sersic profile, with a constant V — / ~ 1.15 at all radii. This 
colour is that of a 10 Gyr old, Z = 0.008 population with 
a Salpeter IMF and an exponentially declining star forma- 
tion rate r = 1 Gyr (SFR oc ex p(— t/r)) predicted by inte- 
grated the iMarigo et al.l (120081 ) isochrones and reweighting 
the timesteps according to the star formation history. Due 
to the relatively high metallicity assumed, this represents an 
upper limit on the V — I colour that any normal halo pop- 
ulation can be expected to display. The plotted star counts 
profile is described by a Sersic profile: 



Hv ~ lJ-v,o + 



2.5b 
log(lO) Vreft 



l/n 



(1) 



with Sersic index n = 5, fj,v,o = 17 mag arcsec"^, b = 2n — 
0.324 and red — 10.0 kp c, resulting in an J -band profile 
similar to that derived bv lZibetti et al.l (|2004l ) in the polar 
direction away from their stacked edge-on disks. 

The effects of EBL extinction on surface photometry 
(blue and red solid lines in Fig. [2J;) have based on the dif- 
fuse EBL levels derived in Zackriss on ct al . ( 2009i) , a Milky 
Way extinction law and a radially constant V-band halo 
opacity of Ay = 0.05 mag. The corresponding surface pho- 
tometry prediction in the case of ERE (blue and red solid 
lines in Fig. [2]d) has been derived assuming that the dif- 
fuse ERE flux is limited to the J-band, while requiring the 
same colour gradient (shown in Fig. [2]i) as the EBL extinc- 
tion scenario. The bottom-heavy IMF scenario prediction 
is based on the assumption that the standard halo at large 
radii is gradually overtaken by a second halo population with 
V — I = 1.5 , as would be expected for a. Z — 0.008, 10 Gyr 
old r = 1 Gyr population with an IMF slope a = 4.50 
(diV/dM a M~" throughout the mass range 0.08-120 



one of the mo dels discussed by [B crgvall et al.l (|201C 
Based on the same lMarigo et al ] (|2008l ) isochrones as we use 
for the standard halo, we predict that the difference between 
the /^totai — ^counts offsets of the standard and bottom- heavy 
IMF halos are A^v —0.75 and A^i « —1.1 in this case. 
The scaling of the relative fluxes of two components are then 
determined by the requirement of producing the same inte- 
grated V-I gradient (Fig. [2ji) as the EBL extinction and 
ERE scenarios. 

While these choices of model parameters and profiles 



may seem arbitrary, we stress that our strategy for pinning 
down the mechanism responsible for the red excess is in- 
sensitive to the details of the adopted surface brightness or 
extinction profiles. Instead, the proposed test is based on the 
offset between /^totai and /icounts over the limited isophotal 
range where a) surface photometry has revealed anomalous 
colours, and b) both /Xtotai and /icounts can be measured (see 
Sect 4.1 and 4.2). In this case, it is only the typical val- 
ues of /itotai, /icounts and extinction within this region of the 
halo that matter. The appearance of the profiles outside this 
region are not important. While the well-known structural 
components of galaxies (thin disk, thick disk, bulge and stel- 
lar halo) can be combined in endless ways to produce the 
same surface brightness level at a given radius, the offset be- 
tween /itotai and /icounts at that point still remains a sensitive 
probe of exotic ingredients entering the mix (Fig. [Sj. 

3.3 Choice of filters 

Even though filters other than V and / can in principle 
be considered for the optical part of these measurements, 
these filters appear to be the most suitable, both because 
similar passbands are routinely used in halo observations 
based on star counts and also because surface photome- 
try can be pushed to very faint isophotes at these wave- 
lengths. While the R fllter may seem more attractive than 
/ for doing surface photometry (due to a fainter sky back- 
ground and less f ringing), this st r ategy - used in the halo 
measurements by Jablonka et al. (2010) - would introduce 
additional ambiguities since a red excess detected inV — R 
could be caused by nebular emission (see Sect. [22]) due to 
the location of the Ha emission line in the R fllter. It is 
moreover u nclear if the red exces s detected bv lZibetti et al.l 
(2004^ and iBergvafl et"al] (j2010l ) in r - i really would be 
detectable in V — R (even after correcting for the median 
redshifts of z = 0.05-0.06 of these stacked samples). As far 
as we know, the only paper that lends any support for an 
anoma lous V — c olour in the halo of an edge-on disk is 
that of iRudv et al] (|l997l ) on NGC 5907. However, Rudy et 
al. arrive at their V — 7? proflles by compiling inhomogeneous 
data from the litera ture. While their l/-band proflles (from 
iLeaueux et~ai]|l99dl ) differ o n the tw o sides of the halo, the 
_R-band proflles (from Sacko tt et al.l[l994 ) have been aver- 
aged over both sides of the halo. If data on one side of the 
galaxy is unrel iable, this procedure may result in spurious 
effects. Indeed. [Leaueux et al.l (|l998t ) note that the process- 
ing of the data on the west side of this galaxy (where a bright 
star happens to be located) is more troublesome than on the 
east. They also flnd that the V — R colour they derive by 
combining t heir own BVI measurements with the _R-band 
proflles from lSackett et al] (| 19941 ) produce V — R colours for 
the halo that are seemingly inconsistent with B — V and 
V — I at the same radius. 



4 OBSERVATIONAL REQUIREMENTS 

The need to combine surface photometry data with star 
counts imposes a number of constraints on the observational 
strategy and the selection of suitable targets. 

Optical/near-infrared surface photometry measure- 
ments can easily trace the light distribution in the inner 
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regions of galaxies, but become increasingly difficult as one 
moves outwards, towards fainter isophotes. The depth to 
which measurements of this type can be pushed is ultimately 
limi ted by one's ability to correct fo r instrumental scatter- 
ing (|de Jonell2008l : ISlater et al.ll2009l') and to subtract off the 



sky with sufficient accuracy (jZheng et al ] |l999l : lMihos" et alJ 
I2OO5I : [ Jablonka et al.l [20101 ). Star counts can be extended 
to much fainter isophotes, where the limit is set by con- 
taminants due to foreground stars and background galax- 
ies. However, star count profiles do not extend to the cen- 
tre of the halo because of crowding problems at high sur- 
face brightness levels, as illustrated by the schematic sur- 
face brightness profiles included in Fig. [2] To implement the 
proposed observational strategy, the overlap of the spatial 
regions for which both surface photometry and star counts 
can be carried out needs to be maximized. In the following, 
we only consider edge-on disk galaxies as potential targets, 
since the detection of a genuine stellar halo will be more 
ambiguous for other morphological types and orientations. 



4.1 Surface photometry 

Even for the haloes of disk galaxies, there are substantial 
variations in the reported halo surface brightnes s at which 
anomalously red colours turn up. Using the Ma rigo et al.l 
l|2008h model, we predict that a normal, old (age > 6 Gyr) 
metal-poor {Z < 0.008), dust-free stellar population with a 
normal IMF should display V — I ~ 0.9-1.2. In the case of 
NGC 5907, the most famous case of a red halo (see Sect. [5]), 
the observed c olour profile reaches V — I > 1.2 at ~ 26 
mag arcsec"^ (jLeaueux et al.lll996l . [l998l ). 

In the haloes of stacked high surface bright- 
ness (|Zibetti et al.l l2004l l and low surface brightness 
ijBergvall et al.ll2010l ) disks, the anomalously red r — i colours 

28- 



iu\V~I measurement (error a(y — I) ^ 0.15), the ay and 
ai errors on the surface brightness measurements at these 
isophotes should not exceed 0.1 mag. Optical surface pho- 
tometry on si ngle targets have been pushed to similar lev- 
els before fe.g. [Barton fc Thompsonll997l : IZheng et al.lll999l : 



are reported at fii ~ 26.7 mag arcsec^ and fii 
29 mag arcsec"'^ respectively. Converting the latter i-band 
results into the 7-band requires knowledge of the spectral 
energy distribution of the halo in the « 6500 — 11000 A 
wavelength range, but the conversion term is - at least for 
a halo dominated by direct star light - not very sensitive 
to the stellar population properties adopted. For a metal- 
poor {Z ^ 0.008), old ~ 10 Gyr population with a stan- 
dard IMF, the conversion is i — / ~ 0.5, whereas bottom- 
heavy IMFs of the type required to fit th e gri halo data 
ijZackrisson et al.ll2006l : [Bergvall et a"l]|2010f ) give i-7 « 0.6. 
Adopting i — I ~ 0.5, these i-band isophotal levels con- 
vert into /i/ ~ 26.2 mag arcsec"'^ and /x/ ~ 27.5-28.5 mag 
arcsec"'^ respectively. 



Hence, while the lZibetti et al.] ()2004h ha lo turns red at 
isopho tes similar to those of NGC 5907, the iBergvall et al.l 
1)20101 ) halo does not turn red until it is traced 1.5-2.5 mag 
fainter. This could be related to the different intrinsic prop- 
erties (luminosity and disk surface brightness) of the sample 
galaxies, but at the current time, this remains unclear, since 
all of these data have not been subject to the same tests for 
systematical problems (instrumental scattering effects and 
sky subtraction residuals). 

To cover all possibilities, the surface brightness mea- 
surements should ideally reach /i/ « 28.5 mag arcsec"'^. To 
robustly rule out or confirm a red excess, the F-band mea- 
surements should therefore ideally reach ~ 29.7 mag 
arcsec"^ (corresponding toV — I — 1.2). To get a meaning- 



I Jablonka et aLlbOloD ." albeit not with an error this small. 

If an anomalously red V — I colour is indeed detected in 
the halo, and the test outlined in Section |3] favours bottom- 
heavy IMF or ERE solutions over EBL extinction, it may be- 
come necessary to add near-infrared data as well. Since sur- 
face photometry measurements of halos are largely limited 
by the brightness of the night sky, and since the night sky 
becomes progressively brighter towards longer wavelengths, 
the observational prospects for matching the optical {V and 
7) halo measurements with near-infrared data appear to be 
maximized in the J-band. A normal, old (age ^ 6 Gyr) 
metal-poor {Z ^ 0.008), dust-free stellar pop ulation with a 
normal IMF should display V - J « 1.5-2.0 (|Marigo et all 
|2008| ). whereas a bottom-heavy IMF would generate a.V — J 
redder than this (i.e. V — J > 2.0). To be able to rule out 
or confirm an anomalously red V — J colour, the fiv ~ 29.7 
measurements would then have to be complemented with 
fij ~ 27.7 mag arcsec"'^ data (allowing a J-band detection 
if 1^ — J ^ 2.0). On the other hand, if the red excess turns 
up already at fiy ~ 27.2 mag arcsec"'^ and /ij « 26.0 mag 
arcsec"^ as in the case of NGC 5907, one could get away 
with fij ^ 25.2 mag arcsec"^ data. Once again, the error 
sho uld ideally not exceed ctj ~ 0.1 mag at this isophote. 

iNoeske et aP (|2003l ) have demonstrated that ground- 
based J-band surface photometry can reach J ~ 25 mag 
arcsec"'^ with a 3.5 m telescope. Assuming that shot noise 
dominates the observations, one should be able to push the 
errors down to the required 0.1 mag by scaling up to a 
bigger groundbased telescope, using longer exposure times 
and integrating over a larger halo region. We estimate, that 
to reach this level while maintaining aj ~ 0.1 mag with 
the 77^ lVi<'-7 instrument at the ESO 8.2m Very Large Tele- 
scop^ one would require about 3 hours of observing time 
(see Appendix A for details). 

However, it is unclear if observations at these isophotal 
levels really are photon-limited, since the J band sky is dom- 
inated by OH emission and water vapor absorption, both of 
which vary on both small and large spatial and temporal 
scales. To reliably push the surface photometry faintward of 
fij ^ 25 mag arcsec"'^, one would therefore need to resort to 
space-based observations. Here, the upcoming James Webb 
Space Telescope ( JWST) 0, tentatively scheduled for launch 
in 2018, may be one of the most suitable choices. When inte- 
grating over 0.37 arcmin^, i.e. a 10" strip across the field of 
view (2.2' X 2.2') of the JWST NIRCam instrument, we find 
that a total of 2 hours in sky-chopping mode should be suf- 
ficient to reach the ideal depth of /xj ~ 27.7 mag arcsec"^. 
The details of these estimates can be found in Appendix A. 
In principle, data of this type can also be used to extend the 
star counts into the J-band, since the angular resolution of 
the JWST is expected to be superior to that of 77S'r in this 
wavelength range. 
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4.2 Star counts 

Star counts in the haloes of nearby galaxies can readily be 
extended to surface brightness levels beyond ~ 30 mag 
arcsec"'^. However, at bright levels, substantial crowding 
problems may occur. To implement the proposed test, one 
would need the star counts in V and I to trace the surface 
brightness profile from regions where no red excess is de- 
tected and into regions where it is. Since the red excess may 
turn up at isophotes as bright as fiv ~ 27.2 mag arcsec"^ 
and Hi « 26 mag arcsec"^ (see Sect. I4.l"|l . one would ideally 
like to start the star counts in regions at least two magni- 
tudes brighter than this, i.e. at levels around /iv ~ 25 mag 
arcsec"^ and /ij ~ 24 mag arcsec"^. 

In Fig. m we present the number of resolved stars per 
square arcsecond at hv ~ 25 mag arcsec"'^ as a func- 
tion of distance to the target halo. These estimates are 
based on ratios of resolved stars per solar luminosity of 
Natalia /Lv,Q ~ 3 X 10""' (tMu solid line), 6 x 10""' (medium 
solid) and 1 x 10""^ (thick solid) for detection thresholds 
1, 2 and 3 magnitudes below the TRGB in the 7-band, 
based on a populatio n synthesis model constructed using the 
iMarigo et al] (|2008l 'l isochrones at Z ^ 0.008 and age ^ 6 
Gyr. Also included are two crowding limits, based on the as- 
sumption that crowding becomes severe once the number of 
resolved stars per resolution element exceeds ~ 0.1. The two 
limits correspond to groundbased observations (thin dashed 
line; assuming circular resolution element with diameter 1") 
and observations with the Hubble Space Telescope (thick 
dashed; assuming circular resolution element with diameter 
0.1"). These assumptions result in crowding limits in rough 
correspondence with the isophotal levels where the incom- 
pleten ess b ecomes seve r e in t he star counts of llrwin et al.l 
l|2005h and ISeth et al.l (|2007l ) for groundbased and space- 
based observations, respectively. 

For groundbased observations, the crowding limit at 
fiv ~ 25 mag arcsec" is reached at distances of D > 1.6 
Mpc (i.e. slightly beyond the Local Group), with only a 
weak dependence on the depth of the observations. For Hub- 
ble Space Telescope (HST) observations, the limit sets in at 
D > 16 Mpc. Since there are very few high-inclination disk 
galaxies within 1.6 Mpc from the Milky Way, the HST ap- 
pears to be the best choice for the star counts at the current 
time. 



5 DISCUSSION 

The need to obtain matching star counts and surface pho- 
tometry data requires a nearby target (Sect. 14. 2p . with a 
large angular halo diameter as a result (~ 5-10'). Since the 
the HST currently represents the best option for the star 
counts (with a field of view of 3.3' x 3.3' for the ACS in- 
strument), the proposed test will realistically be limited to 
a small region of the halo. While this means that halo sub- 
structure (in the form of satellite galaxies and tidal streams) 
may well interfere with the proposed measurements, this 
should not affect the outcome of the proposed test, since 
such structures are expected to affect the star counts and 
surface photometry in similar ways. In fact, the modern view 
of a stellar halo is that of large collection of substructures 
that slowly dissolve into the diffuse component, which is 
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Figure 4. The number of resolved stars per arcsec'^ as a function 
of distance for an old, metal-poor stellar population at /^v ~ 25 
mag arcsec"^. The diagonal lines correspond to the number of 
stars per arcsec^ when the star counts are able to resolve stars 
1 mag (thin solid line), 2 mag (medium solid) and 3 mag (thick 
solid) below the TRGB in the /-band. The dashed horizontal lines 
correspond to the limits for where crowding is expected to become 
serious for star counts carried out from the ground (thin dashed) 
and with the HST (thick dashed). The requirement of being able 
to perform star counts at jiy 25 mag arcsec"^ without seri- 
ous crowding problems limits the distance to the target halos to 
less than 1.6 Mpc (logj^g D(Mpc) ^5 0.2) for groundbased obser- 
vations and 16 Mpc (logjg D(Mpc) ^ 1.2) for HST observations 
(assuming data reaching 2 mag below the TRGB). 

crea ted by both accretion events and in-situ star f ormation 
fe.g. lBullock fc Johnstonll2005l : IZolotov et al.ll2009l ). At this 
time, it is unclear whether the red excess originates in the 
diffuse halo component, in halo substructure or in both. 

Most of the massive (l/rot>100 km/s), edge-on (incli- 
nation > 85 degrees) disk galaxies that have reliable dis- 
tance estimates that meet the distance criteria for HST 
star counts {D ^ 16 Mpc) h ave already been targ eted by 
the G HOST^ survey (|de Jong 2007; Radburn-Smi th et all 
l20lj) . which means that the star counts data are already 
available. To implement the proposed test, one therefore 
simply has to obtain matching surface photometry data. Be- 
low, we discuss a number of potential targets for which star 
counts are already available from GHOSTS or elsewhere - 
NGC 5907, NGC 7814, NGC 3109 and NGC 4565. We argue 
that NGC 4565 is likely to be the best target for this endeav- 
our, and present a worked example of the required exposure 
times on telescopes of various sizes. As will be demonstrated, 
a 4-8 m class telescope represents the most practical option 
for the V and 7-band surface photometry. 

The large diameters of these halos at the relevant 
isophotes (« 5-10') imply that the surface photome- 
try measurements either have to be carried out with 
wide-field imagers (e.g. Subaru/Suprime-Canfl, VISTAQ, 

Galaxy Halos, Outer disks, Substructure, Truncations, Star 
Clusters 

* http://www.naoj.org 
^ http://www.vista.ac.uk/ 
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CFHT/MegaCairff]) - or in sky-chopping mode. Both ap- 
proaches are admittedly very challenging at the extremely 
faint isophotes where the red excess is likely to turn up and 
a combination of both may provide the best route in the 
long run. 

While wide-field imagers save observing time, one tends 
to have better overview and control of the bias level and 
flatfielding accuracy in single chip imaging with single 
amplifier readout. Additionally, the single chip/amplifier 
setup avoids or at least significantly reduces the efi'ects 
of filter-dependent large-scale variations of the pixel re- 
sponse, crosstalk and amplifier gain differences. On the 
other hand, due to the necessary proximity of the tar- 
get galaxy, using a single chip would imply that the level 
of the background sky cannot be estimated from the tar- 
get observations. Therefore chopping mode observations are 
necessary in order to obtain separate imaging of the sky. 
Chopping mode is not a traditional observing strategy in 
the optical. However, it is routinely used to reach sur- 
face brightness level s ~ 10000 times below the sky in the 
near-IR r egime fe.g. 'Cairos et al"2003l: iNoeske et al]l2003l : 
IVaduvescu . Richer & McCall 2006), where the temporal sky 
variations are rapid and can easily result in residual low fre- 
quency variations of the background in the sky-subtracted 
frame. The optical sky is in contrast much more stable and 
there is no obvious reason why chopping mode in the optical 
would not be at least as successful as it is in the compar- 
atively more troublesome NIK regime. For this reason, our 
exposure time estimates are based on the assumption that 
the surface photometry observations are carried out in sky- 
chopping mode. 

5.1 NGC5907 

The edge-on disk galaxy NGC 5907, at a TR GB distance 
of 16.8±0.8 Mpc (|Radburn-Smith et al.ll201lf ). by far rep- 
resents the most famous c ase of a claimed 'r ed halo' (e.g . 



Sackctt ct al. 1994; Lcguc ux et al.l 1996: Ru dv et al 



Leaueux et al.ll 19981 : 1 James fc Casalill 19981 : iZepf et ahkoOOl ) 



19971: 



As discussed in Sect. 14. Il the observed V — I profile reaches 
an anomalous value already at fij 26 mag arcsec"^ 
(|Leaueux et al.iri996l . Il99^ ). Interest in this object declined 
with the discovery of what appeared to be the r emnants of a 
disru pted dwarf gala:xy on one side of the halo (|Shang et al.l 
Il998h and suggestions that this ring-like feature, in combi- 
nation with contamination by foregrou nd stars, could hav e 
resulted in the spurious halo detection (Zheng et al.lll999l ). 
To some extent, this is a matter of semantics, since tidal 
features are now adays regarded as ubiquit ous features of 
stellar halos fe.g. iBullock fc Johnstonll2005h . The details of 
how this tidal feature supposedly explains the extreme halo 
colours reported by others remain unclear, especi ally since 
the r ing itself appears to have normal colours tZheng et al.l 
1 19991 ). 

While the colour anomaly of NGC 5907 is well reported, 
its distance makes it a very difficult target for our test. 
To obtain accurate point source photometry 2 mags be- 
low TRGB would require at least 15 HST orbits per point- 
ing. Several HST pointings would probably be required, as 



iMartfnez-Delgado et al.l (|2008l ) have reported a complex sys- 
tem of streams around this galaxy, and halo measurements 
on and off the streams would be required to determine the 
origin of the colour anomaly. The 6 orbit HST observation 
of the GHOSTS survey was barely deep enough to provide a 
TRGB distance to this galaxy and contains too few stars for 
a reliable luminosity profile. It is furthermore located poorly 
on the major axis near the disk truncation at the end of the 
disk and gives no clean halo measurement. 



5.2 NGC 7814 

The e dge-on galaxy NGC 7814, a t a distance of 14.4±0.7 
Mpc (|Radburn-Smith et al.l I2OIII ) has also been argued 
to exhibit an anom alous V — I colour in its outskirts 
jLegueux et al.l [l995l '). This object is extraordinary in the 
sense that the excess seemingly turns up at relatively bright 
isophotes {fij ~ 23-24 mag arcsec"^). No follow-up obser- 
vations have apparently been undertaken. However, due to 
the very prominent bulge in this galaxy, and the fact that 
the colours were not measured in the polar direction but in 
the outer disk region, it is not clear which component of the 
galaxy is responsible for this. 

The prominent bulge and the relatively large distance 
make NGC 7814 a diflicuh target for our test. The HST star 
counts suffer from crowding at about fij ~ 24 mag arcsec"^ 
at this distance, exactly where the colour anomaly becomes 
significant if the Lequcux ct al. ( 1995) measurements are 
correct. Extensive HST imaging data does exist for this 
galaxy with 4 ACS and WFC3 fields along the minor axis 
out to about 14 arcmin radius (~60kpc) and several more 
fields along the major axis. While absolute matching of the 
star counts and integrated light profile in the central region 
of the galaxy may not be possible, the relative profiles could 
definitely be compare d over at least 3 magnitudes in sur- 
face brightness where iLegueux et al.l (|T995) claim a strong 
colour gradient. However, a cursory look at the pipeline- 
reduced Sloan Digital Sky Surve^Q mosaics (data release 8) 
does not give any obvious support for an anomalous or even 
rising g—i colour at the distance where [Leg ueux et alj (|l995l ) 
note a V — / excess. 



5.3 NGC 3109 

NGC 3109, at a distan ce of 1.32 ± 0.06 Mpc 

iHidalgo. Aparicio fc GallartI |2008| ). represents one of 
the few high-inclination disk galaxies located sufficiently 
close to make groundbased star counts a realistic alternative 
for our proposed test. Star counts reaching « 2 mag below 
the TRGB in the V and / bands are already available 
in regio ns located i n the direction perpendicular to the 
disk CHidalgo et al.l l2008h . allowing the surface densities 
of resolved stars to be traced « 10' (« 4 kpc) into the 
halo (out to fiv ~ 30 mag arcsec"^). There is, however, no 
reported colour anomaly in the halo of this galaxy, and its 
relatively low Galactic latitude {b — 23.1 deg) makes surface 
photometry measurements prone to the effects of Galactic 
cirrus and scattered light from foreground stars. 



http://www.cfht.hawaii.edu/Instruments/Imaging/MegaPrime/ ^ www.sdss.org 
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5.4 NGC4565 

NGC4565 is a near-perfect edge-on galaxy that is ideally 
suited for our research. It is at a distance of D = 11.9 ± 0.3 
Mpc as determined from t he TRGB from HST imaging 
l|Radburn-Smith et al.ll2"oill ). which is near enough that the 
RGB population can be comfortably be measured with HST, 
while it has an angular diameter small enough {D25 ~16 
arcmin) that a large fraction of its halo can be observed 
within a single image of a groundbased camera, allowing 
good flatfielding checks. NGC4565 has a high galactic lati- 
tude (6 = 86.4 deg) and therefore Galactic cirrus and scat- 
tered light from foreground stars will have minimal effect 
on the integrated light observations. There are no obvious 
tidal streams in existing deep images, but an Hi tidal bridge 
exists to a small dwarf galaxy to the North. However, no 
anomalous halo colours have been reported for this galaxy 
as of yet. 

Existing HST images obtained with the Advanced Cam- 
era for Surveys {ACS} and Wide Field Camera 3 { WFC3 
cameras as part of the GHOSTS survey with long expo- 
sure times (>7000 s) in F606W (equivalent to V^-band) and 
F814W (/). So far 5 fields have been obtained, two of which 
along the minor axis, and one more observation with the 
ACS and Wide-Field Camera 3 {WFC3) cameras working 
in parallel has been scheduled. The observations are deep 
enough to reach ~ 2.0 mag below the TRGB, which is suffi- 
cient for the proposed test. Starting a few kpc above the disk 
the color-magnitude diagrams are completely dominated by 
RGB and AGB stars with no young stars indicated. Even 
at 30 kpc above the midplane of the galaxy, the estimated 
contamination of unresolved background galaxies is still less 
than 2%. 

To match these HST star counts with V and 7-band sur- 
face photometry data, new ground-ba s ed ob servations would 
be required. Based on the IWu et al.l (|2002l ) igeeo-band sur- 
face brightness profile (similar to the R band), we estimate 
that the surface brightness level of the NGC 4565 halo in 
the most distant CHOSTS field (at « 6' from the disk) 
to be hr ~ 28 mag arcsec"'^. This can be converted into 
corresponding or vaiues, by assuming a shape for 
the observed halo s pectrum in this wavelength range. The 
iMarigo et al.l l|2008l ) model predicts that an old and metal- 
poor stellar population with a standa rd IMF would display 
V — R ~ 0.5 (see also the discussion in lJablonka et al"]|2010l ) 
and R — I 0.5, but if this were the case, there would 
be no red exces s and the t est would not be me aningful. 
IZackrisson et al" I l|2006l ) and iBergvall et al] (|2010l ) present 
models for old stellar populations with bottom-heavy IMFs 
that are able to reproduce the anomalous gri halo colours 
of stacked edge-on disks. These models predict V — R'^ 0.8, 
R — I « 0.9. EBL extinction models capable of fitting the gri 
data can produce normal g — r colour s yet r— i as red as thos e 
predicted by a bottom-heavy IMF (jZackrisson et al.ll2009l ). 
suggesting F — 7? « 0.5, i? — / « 0.9. However, since the 
J?-band surface brightness is fixed, the observations become 
more challenging if V — J? is red (implying a fainter V-band 
profile) rather than blue (brighter l/-band profile). In or- 
der not to produce estimates that are overly optimistic, we 
therefore adopt the more extreme estimates on both colours: 
1/ - « 0.8 and i? - / 0.9. The ^ 28 mag arcsec"^ 
isophote may then correspond to /ly ~ 28.8 mag arcsec"^ 



and fii ~ 27.1 mag arcsec"^. This is somewhat brighter than 
the ideal requirements derived in Sect. 14. H but sufficient to 
un cover a red ha l o if it is similar to the stacked halo detected 
bv lZibetti et al.l (|2004l '). /n k: 26.2 mag arcsec"^. 

By integrating over 10" strips across this CHOSTS field 
(3.3' X 3.3'), we estimate that one can reach the required lev- 
els {fj,v ~ 28.8 mag arcsec"^, fj,i « 27.1 mag arcsec"^) with 
an error of a = 0.1 mag in approximately 2 hours for the 

V band and 4 hours for the /-band with the 8.2m VLT and 
the F0RS2 instrument (see Appendix A for details). For 
comparison, the corresponding times on smaller telescopes 
(with diameters 7? = 4, 2, 1 and 0.5m, respectively) would be 
T = 9, 37, 148 and 592 hours for the V band and T = 17, 67, 
269 and 1076 hours for the / band. This programme would 
seem best suited for telescopes with diameters in the 4-8 m 
range, since smaller telescopes require excessive amounts of 
observing time. There is, however, otherwise nothing that 
precludes the use of smaller teles copes for measurements 
of thi s type . As demonstrated by iMartmez-Delgado et al.l 
(|2008l . I2OI0I ), tidal features in halos have in fact been de- 
tected at similar isophotes using long observing campaigns 
at very small telescopes {D ^ 0.5m). 

Should the test tend to favour an excess of I band flux 
over a deficit of V band flux as the likely explanation for 
the anomalous halo colour, follow-up J-band observations 
may be required to distinguish between ERE and a bottom- 
heavy IMF. As argued in Sect. 14.11 a normal stellar halo 
should display V — J ^ 1.5 — 2.0. This could also be the 

V — J colour expected in the case of ERE, since neither 
the V or J filters are expected to be strongly affected by 
this mechanism. To rule out a normal V — J colour in the 
NGC 4565 field, one may therefore need to reach as deep as 
the the fij = /ly — 2.0 isophote, which in this case corre- 
sponds to fij ~ 26.8 mag arcsec"'^. The JIVST appears to be 
the best option for this. The estimates provided in Sect. 14. II 
and Appendix A indicate that one should be able to reach 
this level in less than two hours. Despite having a relatively 
small field of view (2.2' x 2.2'), the improved sensitivity of 
the JWST/NIRCam instrument compared to existing HST 
cameras will also allow resolved halo stars to be mapped 
across the halo, in just a fraction of the HST observing time 
spent on existing GHOSTS fields. 



6 SUMMARY 

Attempts to measure the integrated optical/near-IR colours 
of the stellar halos of galaxies through surface photometry 
have revealed halo colours that are at odds with the halo 
properties inferred from the study of bright, resolved halo 
stars. Several possible explanations - a bottom-heavy stel- 
lar IMF, diffuse light from the interstellar medium and ex- 
tinction of extragalactic background light - have been pro- 
posed, but a test capable of differentiating between these 
is still lacking. We argue that these different scenarios can 
be tested by combining star counts in the V and I bands 
with similar surface photometry data, since the proposed 
explanations should result in different offsets between these 
profiles. When using the HST for the star counts, this test 
can be applied for galaxies up to D « 16 Mpc. Matching V 
and I surface photometry data can be obtained with existing 
ground-based telescopes. We discuss a number of potential 
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edge-on disk galaxies for which halo star counts are already 
available, and argue that the NGC 4565 is likely to be the 
best target for this endeavour. 
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APPENDIX A: ESTIMATES OF REQUIRED 
EXPOSURE TIMES 

The need to target galaxies at distances < 16 Mpc implies 
that the halos of large galaxies will measure many arcmin- 
utes, possibly even degrees, across (e.g. a halo with radius 50 
kpc at a distance of 16 Mpc will have an angular diameter 
of 21.5'). While having a telescope field of view sufficiently 
large to cover the entire halo would certainly be an asset, 
the proposed test in principle requires only surface photom- 
etry observations of the part of the halo in which the star 
counts have been carried out. In the following, we assume 
that the star counts will be space-based, and adopt an area 
of 3.3' X 3.3' (corresponding to the HST/ACS observations). 
Since wide-field imagers tend to be rare on the largest tele- 
scopes, we therefore base our exposure time estimates on 
the premise that the surface photometry needs to cover this 
area only. 

The observational and reduction strategies will there- 
fore involve chopping mode and consequent pair-subtraction 
of images. Since this is not standard practice for optical ob- 
servations, optical exposure time calculators (ETCs) do not 
(unlike NIR ETCs) account for the increased noise that re- 
sults from the pair subtraction. This implies that any expo- 



sure time estimated for on-source observations in the optical 
must be quadrupled to obtain the total required telescope 
time. For the duration of this paper all presented optical 
exposure times are total required times which include on- 
source and off-source observations and a further factor of 2 
to compensate for the pair-subtraction. For HAWK-I, a NIR 
instrument, the total increase factor of the ETC estimated 
time is only 2. 

Al Optical surface photometry 

In the optical, ground-based telescopes should suffice 
since the atmosphere is almost transparent in this wave- 
length range. The estimated surface brightness level of the 
NGC 4565 halo at the distance of the GHOSTS field at « 6' 
or « 20 kpc away from the centre is nv « 28.8 mag arcsec"^ 
and « 27.1 mag arcsec"^. For each radial data point 
along the surface brightness profile, the integration area is 
assumed to have length 3.3 arcmin parallel to the plane of 
the disk and width of 10 arcsec in the perpendicular direc- 
tion. Using the F0RS2 instrument at the 8.2 m VLT, this 
area would be covered by 3.2 x 10* pixels. The total error 
budget here consists of flatfielding errors on the order of 3% 
of the background (cff), low- frequency spatial variations 
of the background across one frame {(Tsky), temporal vari- 
ations in the background from frame to frame due to OH 
emissions from the upper atmosphere (ctoh) and shot noise 
due to the random statistical fluctuations of the photon flux 
arriving at each pixel {asDOAi)- The latter error term is sim- 
ply the standard deviation of the mean (SDOM) of the flux 
over an integration bin. apF is a composite error, which 
includes residual background unevenness after flatfielding, 
shutter shading effects and the effects of sensitivity changes 
as a function of colour across the CCD. 

The OH emissions from the ionosphere do not con- 
tribute significantly to the V band flux, therefore the 
total error in this fllter is assumed to be av = 
-s/'^FF + '^iky + ""sDOM where we set an upper constraint 
on the total error by treating asky and osdom as indepen- 
dent of the flatfielding error. As we approach the near-IR 
range, temporal variations of the OH emission will cause a 
non-negligible increase of the total error budget. This effect 
is already visible in the optical / band and we account for it 
by including an extra term in the total error for that filter: 
(Ti = (T%p + al^y + (Jqjj -I- crgj^ioM where aoH reflects the 
change in the background flux due to the variation of the 
OH emissions from frame to frame. To estimate the size of 
this term we obtained a temporal sequence of surface bright- 
ness variations of the background in the / band by observing 
a relatively empty sky field for the duration of 90 mirmtes 
on a clear photometric night with the 2.52 m Nordic Op- 
tical Telescope and the ALFOSC instrument. Although the 
baseline is not particularly long, a distinct periodic shape of 
the variations was detected. A typical variation of the OH 
emission in one period was measured to be 2% of the back- 
ground - a factor 24 greater than the shot noise in the same 
data set. The aoH term can in principle be kept small if one 
has access to an instrument with very fast readout times, 
but for our immediate purpose we assume the term cannot 
be neglected. Allowing a maximum total error per filter of 
0.1 mag we arrive at a required S/N per pixel of 0.076 in 
the V band and 0.175 in the / band. To reach the surface 
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brightness levels of jj.v ~ 28.8 mag arcsec and fii ~ 27.1 
mag arcsec"^ discussed in Sect. [5] we would need 2.2 h for 
the V band and 4.0 h for the I band with VLT/F0RS2. 

A2 Near-IR surface photometry 

If the red excess is assumed to appear at isophotes as bright 
as fj,j ^ 25.2 mag arcsec"'^, the required observations can 
be performed from the ground assuming they are photon- 
limited. Groundbased observations in the NIR regime are 
completely background-dominated, so the readout noise can 
be safely ignored. To comply with the requirement that the 
total error in each filter should not exceed 0.1 mag over each 
integration area of 3.3 arcmin times 10 arcsec, we thus need 
a signal-to-noise (S/N) per pixel of '--^ 4 x 10~*. With the 
8.2 m VLT and the HAWK-I instrument, we can then reach 
fij ~ 25.2 mag arcsec"'^ in a total of ~ 3 hours. To push 
the surface brightness to even fainter levels we would, how- 
ever, need to resort to space-based observations. The 6.5 m 
JWST, tentatively scheduled for launch in 2018, appears to 
be the most efficient choice. To estimate the exposure times 
for the NIRCam instrument (field of view 2.2' x 2.2'), we es- 
timate the increase in S/N due to the integration area (2.2' 
times 10") and compare to the available S/N measurements 
per resolution element for point sources published on the 
NIRCam website. We assumed a total system transmission 
of ~ 0.4, a readout noise of 15 electrons per readout, and 
a maximum of 20 readouts (destructive, array is reset) for 
the chopping mode. The readout noise per readout clearly 
dominates the total background noise over the integration 
bin as estimated from the available wavelength-dependent 
scattered light measurements at the NIRCam website, so 
we have neglected the background noise contribution to the 
total error budget (the readout noise for one readout is a 
factor of ~ 50 greater than the background scattered light 
noise in the J band). The total error we obtain is therefore 
much smaller than the required 0.1 mag per filter, so we take 
S/N 10 to be conservative and arrive at a maximum total 
exposure time of 2 hours in order to reach surface brightness 
levels of 27.7 mag arcsec"^ (Sect. BT]) in the J band. 



